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1. Introduction 

Mixed action simulations became popular in recent years as they combine the simulation ad- 
vantages of a simple sea quark action with the exact or near exact chiral symmetry of overlap or 
domain wall valence quarks. The price to pay, in addition to an internal inconsistency (unitarity 
violation), is the complication in the analysis. One option is to derive and use partially quenched 
mixed action chiral perturbative formulae. Alternatively, one can match the parameters of the va- 
lence and sea quark actions as well as possible and deal with any remaining difference as part of 
the lattice artifacts. This approach is useful if the chiral perturbative formulae do not exist or the 
numerical data does not allow the fitting of all the parameters, or if one desires more insight into 
the physics contained in a particular set of gauge configurations. 

The effectiveness of the latter approach was illustrated in Refs. ^ ^, where we showed that, 
at least within the 2-dimensional Schwinger model, mixed action simulations with overlap valence 
quarks on rooted or unrooted staggered sea quark configurations reproduce the full dynamical 
overlap results if the overlap mass is tuned appropriately. In this work we report our first results 
along the same lines in 4-dimensional 2-flavor QCD. We show that configurations generated with 
2-flavor staggered quarks at a single lattice spacing but at four different quark masses are consistent 
with 2-flavor QCD configurations, at least for the three heavier masses. We also show how the 
topological charge distributions can be used to determine the best overlap valence matching masses 
and that with these mass values the data sets predict a consistent value for the chiral condensate of 
2-flavor QCD. All details can be found in Ref. 

2. Strategy and Simulation setup 

Our sea quark action is the 2-flavor Asqtad staggered action [@-|6|]. We have generated four 
configurations sets, each consisting of 400-500 12^ lattices at a lattice spacing of about a = 0.13fm. 
The details of the sets are summarized in Table |l|. The level of taste breaking, the ratios of the 
heaviest and hghtest pion masses, is approximated from corresponding 2-1-1 flavor results The 
last column lists the separation of the configurations in terms of unit length molecular dynamics 
trajectories. 

Our valence action is an overlap action based on an improved Wilson kernel on HYP smeared 
links. This action was used in recent overlap simulations ^. 

To investigate if the rooted staggered gauge ensembles are consistent with two-flavor con- 
tinuum QCD one has to consider observables that are sensitive to the vacuum and do not depend 
strongly on the valence quark mass. Spectral quantities are not appropriate, but the low lying in- 
frared eigenmodes of the massless valence Dirac operator offer a good choice. Another quantity 
which we will consider is the topological charge of the configurations. On each configurations we 
ask what valence quark mass matches the staggered configurations the best, i.e. what valence mass 
minimizes the lattice artifacts, the difference between lattice data and continuum QCD. 

3. Eigenvalues of the Dirac Operator and Random Matrix Theory 

Random matrix theory (RMT) captures the universal chiral properties of QCD and predicts the 
distribution of the physical (infrared) eigenvalues of the massless Dirac operator in the £-regime. 
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Taste breaking 


Time separation 
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7.18 


0.01 


3.84(6) 


60% 


5 


M 


7.20 


0.02 


3.82(3) 


34% 
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H 


1.22 


0.03 


3.60(4) 


24% 


10 


E 


7.24 


0.04 


3.64(3) 


18% 


15 



Table 1: Parameters of the nf — 2 staggered background configurations. The molecular dynamics time 
separations between the configurations reflect the autocorrelation of the topological charge. 



The predictions are given in fixed topological sector v and depend on the low energy constant Z, the 
infinite volume chiral condensate. The distribution of the (microscopically rescaled) eigenmode 
XTV is given as 

/'v,«(A£y)=Av,„(mIV;?2/), (3.1) 

where m is the quark mass of the configurations (sea quark mass), which in our case is an overlap 
quark mass that corresponds to the background configurations that were generated by staggered 
quarks, i.e. m is the matching quark mass as described in Sect. ^. The value of m is not known a 
priori and therefore Ay^n depends on two variables, M = niLV and L. We fit the measured eigenvalue 
distribution to random matrix theory at fixed M and predict the chiral condensate £. The systematic 



deviation of the data from the RMT prediction of Eq.(|3.1|) characterizes the lattice artifacts, both 
from discretization errors and from the non-locality of the action. This deviation is the measure of 
consistency between the lattice action and continuum QCD and replaces the residue used in Ref.[|l]] 
for the same purpose. If the rooting procedure is correct, it should scale to zero as the continuum 
limit is approached at fixed physical (matching) quark mass, assuming the simulations are done in 
the region where the RMT predictions are valid. 

To fit the cumulative distributions, we use the Kolmogorov-Smirnov (KS) test that minimizes 
^max' maximal deviation between the measured and the predicted cumulative distributions [O- 



130 . An advantage of the KS test is that there is an explicit and simple form for the confidence level 
of the fit. For a given sample length this quality factor 2ks is a monotonically decreasing function 
of Dmax that gives the probability that the measured distribution is consistent with the analytical 
one. The KS fit maximizes the quality factor 2ks or the product of quality factors if more than one 
distribution is used. 

However, 2ks will go to zero exponentially with increasing statistics if the measured distri- 
bution is not exactly described by the analytic form. In any lattice calculations there are lattice 
artifacts and finite volume effects, so the analytic form is never exactly reproduced, the quality 
factor vanishes as the numerical statistics increases. In the following we fit our data by maximizing 
the quality factor (or products of quality factors) according to the KS test, but describe the goodness 
of the fit by the value D^ax itself to enable a comparison with other results. 

Fig. |l] shows a typical fit of the cumulative distribution for the M (arrist = 0.02) data set at 
M= 13.5. The left panel corresponds to the v = 0, the right panel to the v = 1 sector. Only the first 
modes of the v = and 1 topological sectors are included in the fit. In addition to the two fitted 
modes we also show the non-fitted second modes in the same topological sectors. D^ax is almost 
a factor of two smaller for the « = 1, v = 1 mode, but not significantly worse for the non-fitted 
modes than for the fitted n = I, v = mode. 
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Figure 1: RMT predictions of the cumulative distribution of the two lowest eigenmodes in the v and 
1 sectors of the M set at M= 13.5 (see below). The fit uses only the first mode in each topological sector. 
Arrows indicate the maximal deviation between the data set and the analytical predictions. 



In Fig.^ we plot the maximal deviations D^ax as a function of the RMT parameter M. Ev- 
idently the quality of the fit is not very sensitive to the parameter M. While small values are 
disfavored, larger values are almost equally probable. Contrary to our original hopes the eigen- 
mode distributions cannot be used to define a matching mass, it defines only a range of acceptable 
values. 

Result of the fit are similar for the other 
three data sets. The upper panels of Figure || 
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show Dmax fof the fitted modes, and the depen- 
dence on the staggered mass is obvious. Dmax 
is significantly lower at the heaviest E data set 
than for the lightest L one, with the interme- 
diate mass sets lying in between. This behav- 
ior is expected since at finite lattice spacing a 
smaller staggered mass leads to increased taste 
symmetry breaking (Table |l]), it differs more 
from the flavor symmetric valence quark sec- 
tor. With decreasing lattice spacing at fixed 
physical quark mass this deviation should de- 
crease and eventually vanish in the continuum 
limit. 

At each M value the fit predicts YV /a and using ro/a from Table [l] this can be converted to 
physical units as shown on the lower panels of Fig.^ The corresponding overlap mass values m 
are shown along the upper border of the figure. In order to predict the chiral condensate we have 
to find an independent quantity that predicts the matching valence quark mass. The topological 
charge distribution is a possible choice as we will discuss in Sect. ^ 

In Table ^ we list the number of configurations, the D^ax values of the RMT fit at specific 
M = mLV values and the corresponding quality factors. The D^ax values can be compared to those 



Figure 2: Dmax as the function of M for the M data 
set. The fit uses the first modes of the v = and 1 
sectors (filled points). 



from Ref. [13]. That work uses dynamical overlap configurations at similar physical volumes at 
slightly coarser lattice spacing. Using the same fitting strategy as ours they find D^ax = 0. 1 1 ~0.20. 
In view of these numbers we can conclude that, as far as the Dirac operator eigenmode distribution 
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Figure 3: Dmax and E'/^ in GeV as the function of M = niLV for all four data sets. 



is concerned, the rooted staggered action configurations do not show larger lattice artifacts than the 
overlap ones. 



4. Topology 

Since we have sufficient statistics, over 400 approximately independent configurations at each 
coupling value on not too large volumes (12^ or about 6 fm"*), we can study the topological charge 
distribution. Following the discussion of Refs. [|l^-^], we write the probability of encountering 
a charge itv configuration in the dynamical ensemble as Py = Zy{m'LV)Qy{a). Here Qy is the 
quenched probability of a charge ibv configuration, expected to be Gaussian up to 1 /V corrections, 
while Zy describes the suppression due to the fermionic determinant. The fermionic suppression 
factor has been calculated both within chiral perturbation theory and the random matrix model 
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85 
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193 
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Table 2: Results of the RMT fit to the lowest eigenmodes in the v = 0, 1 sectors. For the determination of 
M see Sect. ^. The first error on the matching overlap mass am is due to the uncertainty of E only, while the 
second one takes into account both the errors of E and M. 
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[14, |T^]. Thus the charge probability distribution Py depends on two variables, M = mLV and o. 
The latter can be determined from the quenched topological susceptibility, so a one parameter fit to 
the topological charge distribution data predicts M, which we list in Table 

5. The chiral condensate and matching masses 

With the M = miy values predicted from the topological charge distribution we are now able 
to extract the physical value of the chiral condensate. Combining the M values with r^/a from 
Table |] we find that all four configuration sets predict a consistent value for the £ condensate, as 
listed in Table ^. The only sign that the light L set differs from the RMT prediction more than the 
other mass values is the larger error of the predicted condensate. The value we obtain, 

e;/^^= 291(5) MeV, (5.1) 

is the lattice condensate. It is consistent with predictions obtained on overlap dynamical configu- 



rations []13|], further supporting our observation that the rooted staggered configurations are QCD 
like, the non-local terms of the action can be simply taken into account as lattice artifacts. 

To connect the value of the condensate to a more conventional scheme, like MS at 2GeV, 
one needs the corresponding renormalization factor Z^. Such a factor should be calculated non- 
perturbatively on the staggered configurations with our specific valence Dirac operator. We have 



not done this calculation yet but similar ones exist [12, 13, 17 1. seems to be largely independent 



of the detailed properties of the background configurations and we estimate its value to be Z^ > 0.9, 

1/3 

which will lower Zj^^ by 3%. In addition, there is a finite volume correction to the condensate that 



could lower its value further []13|]. These effects will have to be investigated but they are beyond 
the scope of the present work. 

Combining M and £ we get the values listed in the "am" column of Table ^. These matching 
masses are not only surprisingly large but they do not depend linearly on the staggered masses. 
While the staggered quark mass changes a factor of four between the lightest and heaviest data 
sets, the matching overlap masses change only 50% . This is similar to what we observed in the 
Schwinger model [|l]]. The matching valence masses show an overall shift compared to the stag- 
gered sea mass values. In addition at very small sea quark masses, where the matching breaks 
down, the valence quark masses are largely independent of the sea quark mass values. This is il- 
lustrated in Fig. 3 of Ref. [[I]]. Such behavior implies that staggered configurations at small quark 
masses are not necessarily closer to chiral continuum QCD than the heavier mass configurations. 
All the computational efforts creating light configurations might be in vain, creating only configu- 
rations with larger lattice artifacts. This might not be a problem when the data is analyzed with the 
whole machinery of staggered partially quenched chiral perturbation theory but should be consid- 
ered when individual configuration sets are analyzed in mixed action simulations. Of course this is 
only a lattice artifact and any such effect will disappear as the continuum limit is approached. 

6. Conclusion 

We have studied the properties of the rooted staggered action in a mixed action simulation 
using overlap valence quarks. By comparing physical quantities that are independent of the va- 
lence quark mass to continuum QCD predictions we identified lattice artifacts and studied their 
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dependence on the sea quark masses. In this work we considered the eigenvalue distribution of the 
massless Dirac operator and the distribution of the topological charge. We compared the former 
to the universal predictions of random matrix theory and found that the systematic deviation of the 
data from the predictions were comparable to dynamical overlap simulations. Using the topologi- 
cal charge distribution we could identify the matching overlap valence quark mass value which best 
describes the staggered configurations. We found these matching values to be fairly large and their 
dependence on the staggered sea mass values indicate a finite offset, in addition to a linear mass 
renormalization factor between the valence and sea mass values. With the use of this matching 
mass we extracted the value of the chiral scalar condensate. We found that the predictions from 
all of our staggered configuration sets were consistent. These findings indicate that at our lattice 
spacing, a^iO.lSfm, and with not very light sea quarks the rooted staggered lattice configurations 
have lattice artifacts similar to other lattice actions, the non-local terms arising from the rooting 
procedure can be simply considered as part of the cutoff effects. In order to show that these non- 
local terms indeed become irrelevant in the continuum limit the calculation have to be repeated at 
different lattice spacings and the scaling of the lattice artifacts should be investigated. It would also 
be important to study in a similar manner the lattice artifacts of other observables. 
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